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Abstract 



We give a quantitative account of the ground-state properties of the dilute 
magnetically trapped 87 Rb gas recently cooled and Bose-Einstein condensed 
at nanokelvin-scale temperatures. Using simple scaling arguments, we show 
that at large particle number the kinetic energy is a small perturbation, and 
find a spatial structure of the cloud of atoms and its momentum distribution 
dependent in an essential way on particle interactions. We also estimate the 
superfluid coherence length and the critical angular velocity at which vortex 
lines become energetically favorable. 
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In a remarkable experiment, Anderson et al. |l] have cooled magnetically trapped 87 Rb 
gas to nanokelvin-range temperatures, and observed a rapid narrowing of the velocity distri- 
bution and density profile, which is interpreted as the onset of Bose-Einstein condensation. 



||, in which particles obeying Bose-Einstein statistics condense at low temperatures [|]||; 
indeed, the condensation of trapped atoms has been a "Holy Grail" of atomic physics ||. 
In this paper we give simple quantitative arguments, taking into account effects of the re- 
pulsive interatomic interactions, to determine the properties of the quantum ground state of 
the trapped 87 Rb system, including its geometry, momentum distribution, coherence length 
and critical angular velocity for vortex formation. 

In the experiment, the gas is magnetically trapped in an effective three-dimensional 
harmonic well cylindrically symmetric about the z-axis, with tunable angular frequencies 
u;° in the axial (z) direction and ujj_ = lu®/\/8 in the transverse (xy) plane. The oscillators 
are characterized by lengths a± = (/z/mu/[) 1//2 and a z = (h /mcu®) 1 / 2 , where m is the atomic 
mass. During the condensation phase, with uj° z /2ti = 208 Hz, and a± « 1.25 x 10~ 4 cm, 
the distribution rapidly sharpens with falling temperature, as a macroscopic number, Nq, of 
the Rb atoms begin to occupy the lowest mode of the well. In the absence of interparticle 
interactions the lowest single-particle state has the familiar wave function, 



where f± is the component of r in the xy-plane. The density distribution, p(f ) = N <j) (r) 2 , 
is Gaussian. However, interatomic interactions strongly modify the particle structure in the 
well. 

The low-energy interactions between polarized 87 Rb atoms are repulsive, and are de- 
scribed by an s-wave triplet-spin scattering length, a, determined to be in the range 
85ao < a < 140ao, where ao is the Bohr radius [/J. In the limit in which the density 
varies slowly on a scale a, the interaction energy of the gas per unit volume is given by 
Ei n t = (2Tch 2 a/m)\p(f)\ 2 . The repulsive interactions favor a reduction of the density from 



Trapped atom clouds are new systems, beyond liquid 



4 He and excitons in semiconductors 
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the free particle situation, A^0o(^) 2 - As the number of particles increases, the first effect of 
interactions is to cause the cloud of particles to expand in the transverse direction, where 
the restoring forces are weaker. With further increase in the number, the cloud expands in 
the z-direction. The eventual size of the cloud is determined, in the limit in which the inter- 
particle interactions dominate, by a balance between the harmonic oscillator and interaction 
energies. 

To see the physics of the balance, let us neglect the anisotropy of the oscillator potential 
and assume that the cloud occupies a region of radius ~ R, so that p ~ N/ R 3 ; then the scale 
of the harmonic oscillator energy per particle is ~ vtllo\_R 2 /2, while each particle experiences 
an interaction energy with the other particles ~ (Anh 2 a j m) (N / R?) . The characteristic 
length scale is thus ~ a±(, where the dimensionless parameter characterizing the system is 

C = (87rNa/a±) 1/5 « 4.21 [(a/100a )(A7l0 4 )(l(r 4 cHi/a ± )] 1/5 ; (2) 

under the conditions of the experimental trap with large N, £ ^> 1. The kinetic energy, on 
the other hand, is of order Nh 2 j \2mR 2 ) , so that the ratio of the kinetic to interaction or 
oscillator energies is of order C" 4 ~ iV~ 4/5 ; for N ~ 2000, the kinetic energy in the condensed 
phase is a factor ~ 200 smaller than the interaction energies. 

To estimate the interaction effects more quantitatively we examine the ground state of 
the system in terms of its order parameter ip(r), where / d 3 r\ip(r)\ 2 = N. [We do not 
distinguish N and N Q in the weakly interacting system at zero temperature.] In the Hartree 
approximation, in which ip{f)/N 1 / 2 is the lowest single particle mode, the ground state 
energy of the system in given by a Ginzburg-Pitaevskii-Gross energy functional ||, 



E(ip) = / d 6 r 



-|V^(r)| 2 + i (K) 2 ri + K°)V) W )\ 2 + — |^(r-)| 4 



(3) 
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This approach is familiar in prior studies of Bose-condensed polarized atomic hydrogen 



JT0[1 ; see also Refs. pTT|-[T3[] . The Rb experiments, with lower density, larger atomic mass 



and stronger interactions, fall, however, in a rather different parameter range. 

For a first solution we take ip in the form of the ground-state wave-function, Eq. ([]]): 



with effective frequencies, u± and u z , treated as variational parameters. Substitution of (|J) 
into d||) yields the ground state energy 

* Art (wjl K) 2 u z {u° z ) 2 Nam 1 / 2 1/2 \ 

minimizing E with respect to u± we derive 

co ± = u°JA, (6) 



where 



C 5 (Uz 
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Interactions, by reducing the effective transverse oscillator frequency by A, spread out the 
distribution in the transverse direction by a factor A 1 / 2 ; when iV is sufficiently large that 
C > 1, A 1 / 2 « 2.55 [(iV/10 4 )(a/100a )(10 4 cm/a ± )] 1/4 (w t /u° ± )V*. 

Spreading in the z-direction begins to become significant when the interaction energy 
per particle becomes comparable with hu z ; using (||) and minimizing the resultant ground 
state energy, E{u z ) = Nh (u^A + oj z /A + (o;°) 2 /(4u; 2 )), with respect to u z , we see that this 
condition is (Na/a±)^ J (co z /ujj_) 1 / 2 , which is realized under the experimental conditions. In 
the limit £ 3> 1, the kinetic energy terms in (|5]) are negligible, and we find the shift in the 
frequency in the z-direction, 

2(vrA) 3 / 5 

where A = u®/ujj_] the leading contribution to the energy per particle is 

£ = ^xV^o oc iv 2 / 5 . (9) 
iV 8vr 3 / 5 x v 7 

To obtain the ground state wave function more precisely, we minimize the total energy (|3|) 
with respect to if), keeping the total number of particles fixed, and thus derive the non-linear 
Schrodinger equation 



ft ^2 1 
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where \i is the chemical potential. The physical scales are conveniently brought out by 
rescaling the lengths, letting r = a±^fi, and writing ip(f) = (A^/C 3 a^_) 1 / 2 /(r 1 ), where 
/ d 3 ri\f\ 2 = 1; then ( |10D becomes the dimensionless equation, 

1 



-^ + r^ + A^ + i/(fi)| 



(ii) 



where z/ 2 = (2(i/ C, 2 hu L ). 

In the limit of large N, we can obtain an essentially exact expression for the ground 
state wave function, corresponding to the Thomas-Fermi approximation, by neglecting the 
kinetic energy term, which falls as £ -4 ; then 



f(ri) 2 = v 2 ~ r 2 lx - (X Zl f 



(12) 



in the region where the right side is positive, and / =0 outside this region. This form for the 
wave function is good, except where the density is small, in which case the kinetic energy 
causes the wave function to vanish smoothly. The normalization condition on / implies that 
v = (15A/87I") 1 / 5 , which translates into the relation between fj, and N: 



hu± ( 15A\ 



2/5 



hu ± ( 15AiVa\ 
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(13) 
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Since /x = dE/dN, the energy per particle is simply E/N = (5/7) fi, a result smaller than the 
effective oscillator frequency calculation (|9]) by a factor (36007I") 1 / 5 jl « 0.92. The central 
density of the blob is p(0) = mfx/4:7rh 2 a. 

In the limit of large N, the transverse radius of the cloud is given by R/a± = 
(15AiVa/a ± ) 1/5 , and the half-height in the z-direction is Z = XR. For N = 2000, a = 100a , 
and ujj_/2ir = 208/y/E Hz, one has R/a± ~ 3.24. When the oscillator "spring constants" 
are relaxed in the trap by a factor 75, the equilibrium configuration of the cloud then has 
R/a± m 2.61. These increases are consistent with observations [p]] made when the cloud is 



finally released from the trap | 14}| . We also note that for large N the aspect ratio R/Z equals 



A, whereas in the absence of interactions it is A 1 / 2 ; thus for the experimental conditions, one 
would expect the aspect ratio to be v^8- 

For the Thomas-Fermi wave function, the momentum distribution is proportional to 
| J 2 (k)/k 2 | 2 , where J2 is the Bessel function of order 2, k 2 = (C a _i_) 2 (p± + (Pz/^) 2 ), and p is 
the particle momentum. When Na/a± » 1, the width of the momentum distribution is 
much less than that for a single particle in the oscillator potential. 

The sound velocity, c s , in the interior of the cloud is given by c 2 = (p/m)(dp/dp) = p/m, 
which in the large N limit equals (hui D (15 \N a/ a_i_) 2//5 . In this limit, the lowest mode of 
excitation in the transverse direction of the system has frequency of order R/c s ~ 

The superffuid coherence length, £, which determines the distance over which the con- 
densate wave function can heal, can be estimated by equating the kinetic energy term in 
Eq. (|10D, ~ h 2 /(2m£ 2 ), to the interaction energy, which yields 

e = (STTpa)" 1 , (14) 

where p is the local density. With the central density of the cloud computed in the Thomas- 
Fermi approximation, p(0) = mpj (47ch 2 a), we have 

( f^V i «x y/* (15) 
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Thus when the number of particles is sufficiently large that the Thomas- Fermi approximation 
is valid, the coherence length is small compared with the size of the blob, and the system 
should exhibit superfluid properties more like those of a bulk superffuid than an atomic 
nucleus, where £ ~ R. 

An experimentally important confirmation of Bose-Einstein condensation would be the 
observation of formation of a vortex line in a rotating system. The critical angular frequency, 
fl c i, at which it becomes energetically favorable for a vortex line to be created under rotation 
about the z-axis is 

fi ci ~ -\ 2 ln{ R /&- ( 16 ) 
mR z 



For cloud radii ~ 5x 10 4 cm, this value corresponds to a rotation frequency of order 10 Hz. 
Finally we consider the case of atoms, such as spin-polarized 85 Rb |7| or 7 Li p~T 



with a negative scattering length, corresponding to a low energy attractive interaction. 
A uniform state of such atoms at low density would be unstable to formation of long- 
wavelength density waves, signalling a gas-liquid phase transition. However, as discussed 
theoretically in [|T3?II3, and seen experimentally in ||, the physics in a trap is different; 
this can be understood in the present context by considering the variational calculation 
above. With increasing N, the spatial extent of the wave function is reduced. Provided 
A 2 = (l — (2/ / K) 1 / 2 (N\a\/a±)(u z /u)j_ ) 1 / 2 ^) (cf. (0)) remains positive, the kinetic energy term 
is able to stabilize the system. However, if A 2 becomes negative, the attractive forces over- 
whelm the kinetic energy, and the cloud becomes unstable to collapse. The critical number 
of particles for collapse is ~ (n/2\) 1 ^ 2 a±/\a\. In 85 Ru, for which — lOOOao < a < — 60ao 
0, under the experimental conditions in Ref. [|TJ with uj z /2tc = 208 Hz, this number is 
~ 20 — 300; in the 7 Li trap of Ref. |TT| it is ~ 3000. The final state of the collapsed cloud 
is determined by the shorter-range repulsive components of the interatomic potential. 

To summarize, our calculations provide quantitative results that confirm and extend the 
qualitative considerations in Ref. JTJ on the effect of particle interactions on the properties of 
a cloud of bosons. Experimental confirmation of the dependence on trap parameters, particle 
number and atomic properties of the size of the cloud and the momentum distribution would 
give one increased confidence in the interpretation of the data. 
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